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Summary 

How  dynamic  load  affects  the  pitting  fatigue  life  of  external 
spur  gears  was  predicted  by  using  NASA  computer  program 
TELSGE.  TELSGE  was  modified  to  include  an  improved  gear 
tooth  stiffness  model,  a  stiffness-dynamic  load  iteration 
scheme,  and  a  pitting-fatigue-life  prediction  analysis  for  a  gear 
mesh.  The  analysis  used  the  NASA  gear  life  model  developed 
by  Coy,  methods  of  probability  and  statistics,  and  gear  tooth 
dynamic  loads  to  predict  life.  In  general,  gear  life  predictions 
based  on  dynamic  loads  differed  significantly  from  those  based 
on  static  loads,  with  the  predictions  being  strongly  influenced 
by  the  maximum  dynamic  load  during  contact. 

With  the  modified  TELSGE,  parametric  studies  were 
performed  that  modeled  low-contact- ratio  involute  spur  gears 
over  a  range  of  gear  speeds,  numbers  of  teeth,  gear  sizes, 
diametral  pitches,  pressure  angles,  and  gear  ratios.  Dynamic 
loads  and  pitting  fatigue  lives  were  calculated.  Gear  mesh 
operating  speed  strongly  affected  predicted  dynamic  load  and 
life.  Meshes  operating  at  a  resonant  speed  or  at  one-half  the 
resonant  speed  had  significantly  shorter  lives.  Dynamic  life 
factors  for  gear  surface  pitting  fatigue  were  developed  on  the 
basis  of  the  parametric  studies.  The  effects  of  number  of  teeth, 
gear  size,  diametral  pitch,  pressure  angle,  and  gear  ratio  on 
predicted  life  were  related  to  the  contact  ratio.  In  general, 
meshes  with  higher  contact  ratios  had  higher  dynamic  life 
factors  than  meshes  with  lower  contact  ratios.  A  design  chart 
was  developed  for  use  in  the  absence  of  a  computer  and 
program  TELSGE.  An  example  illustrates  the  use  of  the  design 
chart. 

Introduction 

Gears  may  fail  from  scoring,  tooth  fracture  due  to  bending 
fatigue,  or  surface  pitting  fatigue.  Scoring  failure  is  usually 
lubrication  related  and  can  be  prevented  by  proper  lubrication 
and  proper  operating  temperatures.  Tooth  fractures  are  usually 
caused  by  poor  materials,  improper  design,  or  overloading 
and  can  be  prevented  by  designing  for  bending  stresses  below 
the  material's  maximum  allowable  stress.  The  American  Gear 
Manufacturers  Association  (AGMA)  has  a  standard  practice 
for  predicting  gear  surface  pitting  fatigue  (ref.  1).  The  method 
assumes  that  infinite  life  results  when  the  maximum  surface 
contact  stresses  are  less  than  the  material's  endurance  limit. 
Surface  contact  stress  calculations  may  include  a  dynamic 


factor  to  account  for  gear  dynamic  loading.  The  AGMA 
recommends  a  dynamic  factor  of  1  for  gear  teeth  of  high 
accuracy  but  states  that  actual  dynamic  loads,  computed  or 
measured,  can  be  used  (ref.  1). 

Gear  research  authorities  do  not  completely  agree  on  surface 
pitting  fatigue.  Some  state  that  gear  materials  do  not  have 
surface  endurance  limits  (refs.  2  and  3),  as  is  true  for  rolling- 
element  bearings.  In  1975,  Coy  developed  an  improved  model 
for  the  surface  fatigue  life  of  spur  and  helical  gears,  using  an 
approach  similar  to  that  for  rolling -element  bearings  (refs.  2 
to  6).  This  work  did  not,  however,  include  the  effect  of 
dynamic  load. 

Early  contributions  to  gear  dynamic  loading  were  made  by 
Buckingham,  Tuplin,  Richardson,  and  Attia  (refs.  7  to  10). 
More  recently  computer-based  analytical  programs  have  been 
developed  to  determine  gear  tooth  dynamic  loads  (refs.  1 1  to 
16).  The  dynamic  loads  of  these  programs  depend  on  such 
factors  as  inertia  and  stiffness  of  rotating  members,  tooth 
spacing  and  profile  errors,  size,  and  speed.  The  loads  are 
determined  by  solving  the  equations  of  motion  of  a  given  gear 
mesh  system. 

The  objective  of  the  present  study  was  to  combine  the 
dynamic  load  calculation  procedure  of  Wang  and  Cheng 
(ref.  14)  with  the  NASA  gear  life  model  of  Coy  (refs.  2  to 
6)  to  determine  how  dynamic  load  affects  the  pitting  fatigue 
life  of  external  spur  gears.  NASA  computer  program 
TELSGE,  modified  to  include  Cornell's  gear  tooth  stiffness 
model  (ref.  17),  a  stiffness-dynamic  load  iteration  scheme, 
and  a  pitting  fatigue  life  analysis,  was  used  to  predict  gear 
dynamic  loads  and  life.  Parametric  studies  using  modified 
TELSGE  were  performed  for  low-contact-ratio  involute  gears 
with  no  tooth  spacing  or  profile  errors.  Gear  dynamic  loads 
and  tooth  stiffnesses  were  calculated  as  a  function  of  contact 
position  and  speed.  On  the  basis  of  the  parametric  studies 
dynamic  life  factors  for  gear  surface  pitting  fatigue  were 
developed  as  a  function  of  speed  and  contact  ratio. 

Analysis 

Gear  Life  Model 

Current  theory.  —The  life  model  proposed  by  Lundberg  and 
Palmgren  (refs.  18  to  20)  is  the  commonly  accepted  theory 
for  predicting  the  pitting  fatigue  life  of  rolling-element 
bearings.  Because  the  fatigue  failure  mechanism  is  similar  for 
both  gears  and  rolling-element  bearings  made  from  high- 
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strength  steel,  the  Lundberg-Palmgren  model  for  bearings  has 
been  adapted  to  predict  gear  life  (refs.  2  to  6).  Reference  6 
gives  the  life  for  a  90-percent  probability  of  survival  rj  of  a 
single  tooth  on  a  driver  or  driven  gear  of  a  mesh  as 

n  =  B43/i9Ep-5e-°AQ-43  (1) 

where  B  is  a  material  constant  based  on  experimental  data; 
/is  the  tooth  face  width;  Ep  is  the  curvature  sum  at  the  start 
of  single-tooth  contact;  f  is  the  involute  surface  length  during 
single-tooth  contact;  and  Q  is  the  static  tooth  load,  normal  to 
the  contact.  A  complete  list  of  symbols  is  given  in  appendix  A. 

The  life  of  the  complete  driver  gear  (all  teeth)  Lx  in  terms 
of  driver  gear  rotations  is 

L ,  =  (2) 

where  N j  is  the  number  of  teeth  on  the  driver  gear  and  e ,  the 
Weibull  exponent,  is  a  measure  of  scatter  in  fatigue  life. 
Experimental  research  on  AISI  9310  steel  spur  gears  has 
shown  gear  fatigue  to  follow  the  Weibull  failure  distribution 
with  e  ^  2.5  (ref.  3). 

The  life  of  the  complete  driven  gear  in  terms  of  driver 
gear  rotations  is 

^  (5> 

where  N2  is  the  number  of  teeth  on  the  driven  gear.  The  mesh 
life  (both  driver  and  driven  gears)  L ^  in  terms  of  driver  gear 
rotations  is  given  by 

Lm=  (L,-‘  +  LT'YU'  (4) 

Expanded  theory.— To  adapt  the  current  gear  life  model  for 
predictions  based  on  gear  tooth  dynamic  loads,  the  tooth  was 
divided  into  intervals  (fig.  I).  The  use  of  intervals  allowed 
the  current  gear  life  model  to  account  for  load  and  curvature 
sums  varying  with  contact  position.  The  complete  gear  tooth 
life  was  determined  from  the  interval  lives  and  methods  of 
probability  and  statistics.  The  details  are  as  follows. 

When  a  pair  of  external  spur  gears  is  in  mesh  (fig.  2),  the 
line  tangent  to  the  base  circles  of  both  the  driver  and  driven 
gears  is  called  the  line  of  action.  The  gears  begin  contact  when 
the  outside  radius  of  the  driven  gear  intersects  the  line  of 
action.  As  the  gears  rotate,  the  contact  point  occurs  on  the 
line  of  action.  The  contact  ends  when  the  outside  radius  of 
the  driver  gear  intersects  the  line  of  action.  The  point  at  which 
the  pitch  circles  of  the  driver  and  driven  gears  intersect  is  called 
the  pitch  point.  The  distance  along  the  line  of  action  from  the 
pitch  point  to  the  start  of  contact  is 


Z\  =  r\  i  ~  rp,2  sin  <p  (5) 


Figure  1.— Tooth  intervals  of  a  meshing  gear  pair. 
Driver  gear 


Figure  2.— Basic  geometry  of  a  pair  of  external  spur  gears  in  mesh. 


The  distance  along  the  line  of  action  from  the  pitch  point  to 
the  end  of  contact  is 

Zi  =  VrJ.i  -  rf,  ~  rp  )  sin  *>  (6) 
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The  contact  length  Z  is  defined  as 

Z  =  zx  +  z2  (7) 

Dividing  the  contact  length  into  equal-size  intervals  of  length 
Ax  gives 


(8) 


where  J  is  the  total  number  of  intervals  on  a  tooth,  and 

x i  =  —  Z\  +  (i  —  1 )  Ax  for  i  =  1  to  J  +  l  (9) 

where  jc  is  the  contact  position  along  the  line  of  action.  The 
value  of  x  is  negative  when  contact  is  before  the  pitch  point, 
zero  when  at  the  pitch  point,  and  positive  when  after  the  pitch 
point. 

The  life  of  each  interval  for  a  90-percent  probability  of 
survival  is  given  from  equation  (1)  by 

Vj  =  B4  3/3  sf°  *Qf*  3  for  j  =  1  to  J  (10) 

where  B  and  /do  not  change  from  interval  to  interval.  Both 
curvature  sum  and  involute  length,  however,  change  with 
contact  position. 

At  the  1th  contact  position  the  radii  of  curvature  of  the 
driver  and  driven  gears  (fig.  3)  are 


RU  =  rp.  1  sin  +  x> 

(11) 

*2.i  =  rp,2  sin  <p  -  x, 

(12) 

The  curvature  sum  at  the  Ith  contact  position  is 


&l,i  ^2.i 


(13) 


For  the  7th  interval  the  average  curvature  sum  used  in  the  life 
model  is 


Spy  =  —  +2Sft>l  for  i  =j  (14) 

The  curvature  sum  varied  slightly  with  contact  position  for 
the  example  gear  mesh  data  from  table  I  for  100  intervals  on 
each  tooth  (fig.  4(a)).  (The  contact  position  was  made 
dimensionless  by  dividing  by  the  base  pitch  pb .)  The  plot 
shows  the  curvature  sum  to  be  symmetric  about  the  pitch  point 
(x  =  0).  This  was  true  only  because  the  driver  and  driven  gears 
of  the  example  were  the  same  size. 


Driver  gear 


Figure  3.— Curvatures  of  involute  teeth  in  contact. 


TABLE  I.— BASELINE  DATA  FOR  BOTH  DRIVER  AND  DRIVEN  GEAR 


Number  of  teeth .  36 

Diametral  pitch .  8 

Outside  radius,  cm  (in.) .  .  6.033  (2.375) 

Base  pitch,  cm  (in.) .  0.937  (0.369) 

Face  width,  cm  (in.) .  0.635  (0.250) 

Pressure  angle,  deg . 20 

Root  radius,  cm  (in.) .  5.318  (2.094) 

Fillet  radius,  cm  (in.) .  0.102  (0.040) 

Chordal  tooth  thickness,  cm  (in.) .  0.485  (0.191) 

Normal  load,  N  (lb) .  1718  (386) 

Speed,  rpm .  4000 

Material . Steel 


The  involute  surface  lengths  of  the  driver  and  driven  gears 
for  the  y*  interval  (for  small  Ax)  are 


*2  J 


/Ax \ 

\rb.2/ 


+  Ax  tan  ^  for  f  =  j  (15) 
+  Ax  tan  ^  for  i  =  j  (16) 


The  involute  length  is  a  linear  function  of  contact  position  (fig. 
4(b)).  Equations  (15)  and  (16)  imply  that  rotating  a  gear  mesh 


3 


1  ?5  50  7  5  1  00 

Interval  index,  j 


(a)  Average  curvature  sum  of  interval. 

(b)  Length  of  involute  surface  of  interval. 

(c)  Average  normal  load  on  interval. 

(d)  Life  of  interval  for  driver  gear.  Tooth  life,  13  200  Mrev;  gear  life,  3100 
Mrev;  mesh  life,  2400  Mrev. 

Figure  4.— Effect  of  contact  position  on  gear  life  parameters  for  gear  data 
from  table  I. 

through  equal  angles  produces  unequal  involute  lengths  and 
thus  different-size  tooth  intervals  (as  shown  in  fig,  1). 

By  using  intervals  the  life  model  considers  load  that  can  vary 
with  contact  position.  For  the  7  th  interval  the  average  load 
used  in  the  life  model  is 

Qj  =  -  +2g<*1  for  i  =j  (17) 


The  static  load  variation  with  contact  position  depends  on 
the  number  of  teeth  in  contact  (fig.  4(c)).  As  a  pair  of  teeth 
begin  contact,  the  preceding  pair  of  mating  teeth  are  also  in 
contact.  This  double-tooth-pair  contact  occurs  for  intervals  1 
to  41 ,  and  it  is  assumed  that  half  the  applied  load  is  transferred 
per  contact.  Near  the  pitch  point  single-tooth-pair  contact 
occurs  (intervals  42  to  59),  and  all  the  load  is  transferred  by 
it.  Toward  the  end  of  contact,  double-tooth-pair  contact  again 
occurs  (intervals  60  to  100)  as  the  following  pair  of  mating 
teeth  begin  contact.  As  before,  it  is  assumed  that  half  the  load 
is  transferred  per  contact. 

The  life  of  a  complete  gear  tooth  17,  is  determined  from  the 
interval  lives  and  methods  of  probability  and  statistics  where 

J  \  —  l/e 

E  v')  os) 

)=l  / 

The  complete  tooth  life  was  always  shorter  than  the  lives  of 
the  shortest-lived  intervals  (fig.  4(d)).  Also,  intervals  with 
larger  applied  loads  had  much  more  influence  on  gear  tooth 
life  than  intervals  with  smaller  loads. 

The  tooth  lives  for  a  driver  and  driven  gear  in  mesh  are 
determined  by  the  expanded  life  theory  and  equation  (18).  They 
are  equal  if  the  driver  and  driven  gears  are  the  same  size.  They 
are  slightly  different  if  the  driver  and  driven  gears  are  different 
sizes  because  of  curvature  sums  and  involute  lengths.  The 
complete  gear  lives  and  mesh  life  are  determined,  as  before, 
by  using  equations  (2)  to  (4)  and  substituting  17,  for  77. 

The  total  number  of  tooth  intervals  was  varied  from  30  to 
over  400  to  check  convergence  on  life.  Static  loads  were  used. 
All  cases  predicted  the  same  tooth  life.  Gear  size,  diametral 
pitch,  pressure  angle,  and  gear  ratio  were  also  varied  to 
compare  mesh  lives  predicted  by  the  current  and  expanded 
theories.  Static  loads  were  used.  The  expanded  theory 
predicted  mesh  lives  a  little  longer  than,  but  within  10  percent 
of,  those  predicted  by  the  current  theory  for  meshes  with  equal- 
size  gears.  This  difference  was  caused  by  the  expanded 
theory’s  curvature  sum  variation  with  contact  position.  Thus 
for  meshes  with  equal-size  gears  the  curvature  sum  variation 
had  a  small  effect  on  life.  For  meshes  with  unequal-size  gears, 
however,  there  were  greater  differences  in  the  mesh  lives 
predicted  by  the  two  theories. 

Gear  Tooth  Dynamic  Loads 

Gear  tooth  dynamic  load  model—  The  contact  load  of 
meshing  gear  teeth  varies  as  the  contact  point  moves  along 
the  line  of  action.  This  is  known  as  dynamic  load.  It  is  mainly 
caused  by  single-  and  double-tooth-pair  contact  transitions, 
tooth  stiffness  variation  along  the  contact,  and  tooth  profile 
deviations  from  true  involutes  (tooth  profile  errors).  NASA 
computer  program  TELSGE  (refs.  14  to  16)  was  used  to 
determine  gear  tooth  dynamic  loads.  The  program  models 
meshing  gears  as  a  pair  of  rigid  disks  connected  by  a  spring 
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Driver  gear 


Driven  gear 


Figure  5  —  Dynamic  model  of  meshing  gears. 

(fig.  5).  The  spring  stiffness  corresponds  to  gear  teeth 
stiffnesses. 

The  dynamic  load  model  uses  the  equations  of  motion 
governing  the  angular  displacements  of  the  driver  and  driven 
gears.  By  converting  the  angular  movements  of  the  disks  to 
linear  displacements  along  the  line  of  action,  and  by  algebraic 
manipulation,  the  equations  of  motion  are  represented  by  a 
single  differential  equation,  where 


c, eq  =  2  (23) 

where  f  is  the  damping  ratio,  is  the  equivalent  stiffness 
per  unit  face  width  (discussed  in  the  section  Equivalent  gear 
tooth  stiffness ),  and  Ps  is  the  static  load  per  unit  face  width. 

The  relative  displacement  is  determined  as  a  function  of 
contact  position  by  using  a  Runge-Kutta  numerical  method  and 
solving  equation  (19).  The  dynamic  load  on  a  gear  tooth  is 
determined  as  a  function  of  contact  position  by 

Pd  =  KX  (24) 

where  Pd  is  the  dynamic  load  per  unit  face  width  and  K  is  the 
combined  stiffness  per  unit  face  width  (discussed  in  the  section 
Equivalent  gear  tooth  stiffness ).  Note  that  when  X  is  negative, 
the  teeth  separate  and  the  dynamic  load  is  zero.  Although  tooth 
profile  errors  can  be  accounted  for  in  equations  (19)  and  (24), 
they  were  beyond  the  scope  of  this  study. 

Gear  tooth  stiffness,— Computer  program  TELSGE  was 
modified  to  incorporate  the  gear  tooth  stiffness  model  of 
Cornell  (ref.  17),  regarded  as  the  present  state  of  the  art.  The 
stiffness  model  consists  of  tooth  bending  as  a  cantilever  beam, 
fillet  and  foundation  flexibilities,  and  local  Hertzian 
compression,  all  as  functions  of  contact  position.  In  Cornell's 
model  the  deflections  due  to  bending  and  fillet  and  foundation 
flexibilities  are  expressed  as  linear  functions  of  load,  but  the 
deflections  due  to  Hertzian  effects  are  not  linear  with  load. 
This  makes  the  stiffness  of  a  gear  tooth  dependent  on  dynamic 
load,  and  equation  (19)  nonlinear. 

Equivalent  gear  tooth  stiffness,— The  stiffnesses  of  the 
driver  and  driven  teeth  of  a  mesh,  and  k2  respectively,  are 
found  by  the  methods  of  Cornell  (ref.  17).  The  combined 
stiffness  K  for  a  pair  of  teeth  in  contact  is 


+CatX  +  KeJ(  =  Ps 

(19) 

K  *'*2 

*1  +  *2 

(25) 

The  dependent  variable  X>  called  the  relative  displacement, 
is  the  compression  of  the  spring  along  the  line  of  action. 

For  a  single  pair  of  teeth  in  contact  (fig.  6(a))  the  equivalent 
stiffness  is 

X  —  —  s2 

(20) 

* 

II 

(26) 

where 

For  two  pairs  of  teeth  in  contact  (fig.  6(b))  the  equivalent 
stiffness  is 

5,  =  r*,0,  and  s2  =  rb2$2 

(21) 

'>  +  tfl) 

(27) 

The  equivalent  mass  per  unit  face  width  is 

where 

Af 

*"  Af,  +  Al2 

(22) 

*<■>_  w 

*i'> + k $•> 

(28) 

The  equivalent  damping  coefficient  per  unit  face  width 
includes  the  effect  of  viscous  damping. 

*(n,_ 

(29) 

(a)  Singie-tooth-pair  contact. 

(b)  Double-tooth-pair  contact. 

Figure  6. —Gear  tooth  stiffness  models. 


The  superscript  (I)  refers  to  the  first  pair  of  teeth  in  contact 
and  (II)  refers  to  the  second  pair  of  teeth  in  contact.  The 
equivalent  stiffness  of  equation  (19)  varies  from  double-tooth- 
pair  contact  at  the  start  of  mesh  to  single-tooth-pair  contact 
and  back  to  double. 

Iteration  of  gear  tooth  stiffness  and  dynamic  load.— 
Because  of  the  Hertzian  compression,  gear  tooth  stiffness  is 
not  independent  of  dynamic  load.  TELSGE  was  therefore 
modified  to  iterate  for  dynamic  load  (fig.  7).  First  the  static 
load  is  defined.  As  in  the  example  (fig.  4(c))  all  the  load  is 
transferred  per  contact  during  single-tooth-pair  contact,  and 
half  the  load  is  transferred  per  contact  during  double-tooth- 
pair  contact.  Next  the  combined  stiffness  is  determined  along 
the  contact  position  by  using  the  static  load  in  the  Hertzian 
deflection  computation.  Then  the  dynamic  load  is  determined 
along  the  contact  position.  Next  combined  stiffness  is 
recalculated  by  using  the  calculated  dynamic  load  in  the 
Hertzian  deflection  computation.  Then  dynamic  load  is 
recalculated  by  using  the  latest  stiffness  values.  The  stiffness 
and  load  calculations  continue  until  the  change  in  dynamic  load 
with  each  iteration  becomes  smaller  than  a  preset  amount. 

With  modified  TELSGE  and  the  example  data  (table  I)  the 
dynamic  load  required  only  four  iterations  to  converge  to 
within  0.1  percent  (fig.  8).  So  few  iterations  were  required 
since  the  Hertzian  deflection  was  usually  only  10  to  20  percent 
of  the  total  gear  tooth  deflection.  The  variation  in  equivalent 
stiffness  due  to  double-  and  single-tooth-pair  contact  transitions 
is  a  major  excitation  in  the  dynamic  load  model  (fig.  9).  The 


Figure  7.— Flowchart  of  gear  tooth  combined  stiffness-dynamic  load  inter¬ 
action  scheme  in  computer  program  TELSGE. 


Figure  8.— Effect  of  gear  tooth  combined  stiffness-dynamic  load  interation 
scheme  on  dynamic  load  for  gear  data  from  table  I. 
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tooth-pair  contact 
I _  mean  value, 


<S 


^  Single-  and  double- 
tooth-pair  contact,  Kgg 

/-Single-tooth-pair 
contact  only,  K 


Normalized  contact  position,  x/p^ 


Figure  9.— Effect  of  contact  position  on  gear  tooth  combined  stiffness  for 
gear  data  from  table  I. 


dynamic  load  varied  appreciably  from  the  static  when  the 
operating  conditions  of  the  example  were  used  (fig.  10).  The 
maximum  dynamic  load  during  contact  was  about  30  percent 
greater  than  the  static  load. 
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Figure  10.— Effect  of  contact  position  on  gear  tooth  dynamic  load  for  gear 
data  from  table  I. 


Gear  life  using  dynamic  loads.— The  expanded  gear  life 
model,  which  accounts  for  variations  of  load  and  curvature 
sums  with  respect  to  contact  position,  was  incorporated  in 
modified  TELSGE.  The  dynamic  loads  were  used  in  the  life 
model,  where 

Qi  =  ?d.J  for  i  -  1  to  101  (30) 


(TELSGE  divides  the  contact  length  into  100  intervals.)  For 
the  data  from  table  I  the  mesh  life  based  on  dynamic  loads 
was  then  50  percent  shorter  than  that  based  on  static  loads. 
The  cause  was  the  increase  in  maximum  load  during  contact 
when  dynamic  loads  were  considered  (fig.  10). 


Results  and  Discussion 

NASA  computer  program  TELSGE,  modified  to  include  an 
improved  gear  tooth  stiffness  model,  a  tooth  stiffness-dynamic 
load  iteration  scheme,  and  a  pitting  fatigue  life  prediction 
method,  was  used  to  perform  parametric  studies.  Dynamic 
loads  and  gear  mesh  life  predictions  were  performed  over  a 
range  of  gear  speeds,  numbers  of  teeth,  gear  sizes,  diametral 
pitches,  pressure  angles,  and  gear  ratios. 

Effect  of  Speed  on  Dynamic  Load  and  Life 

Modified  TELSGE  was  run  using  the  mesh  data  in  table  I 
for  speeds  ranging  from  600  to  12  000  rpm.  At  very  low 
speeds  the  dynamic  load  as  a  function  of  contact  position 
(fig.  1 1)  resembled  the  static  load.  However,  spikes  occurred 
at  double-  to  single-tooth-pair  contact  transitions,  and  at  single 


lTxlO3 


1.0  0 


Figure  1 1.— Effect  of  contact  position  and  speed  on  gear  tooth  dynamic  load  for  gear  data  from  table  1. 
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to  double.  As  the  speed  increased,  the  dynamic  load  as  a 
function  of  contact  position  differed  appreciably  from  the 
static. 

The  dynamic  load  reached  a  maximum  at  a  resonant  speed 
o}n  of  about  8500  rpm.  At  speeds  below  resonance  the 
excitation  frequency  from  the  change  in  equivalent  stiffness 
was  lower  than  the  resonant  frequency,  and  the  dynamic  load 
was  basically  an  oscillatory  load  superimposed  on  the  static 
load.  This  produced  peak  dynamic  loads  greater  than  the  static 
load.  At  speeds  above  resonance  the  dynamic  load  had  a 
smoother  response,  with  peaks  lower  than  the  static.  This  was 
caused  by  the  greater  inertia  forces  at  higher  speeds.  The 
resonant  speed  can  be  approximated  by 


cosy  /  60\ 
N  \2tt/ 


(31) 


Here  although  the  mean  equivalent  stiffness  varies  with 
load  and  speed  due  to  Hertzian  effects,  its  influence  on  a>„  is 
not  significant. 

For  the  data  in  figure  1 1  the  maximum  dynamic  load  during 
contact  was  greatest  at  the  resonant  speed  (fig.  12).  It  was  also 
greater  than  the  static  load  at  speeds  below  resonance,  with 
a  secondary  peak  at  about  u)/oo„  =  0.5.  At  speeds  above 
resonance  the  maximum  dynamic  load  during  contact 
decreased  and  was  less  than  the  static  load  above  =  1.2. 

The  gear  mesh  life  as  a  function  of  speed  for  the  mesh  data 
in  table  I  is  shown  in  figure  13.  The  dynamic  life  factor  is 
defined  as 


Cv  =  y  (32) 

where  Ld  is  the  gear  mesh  life  based  on  the  expanded  life 
theory  and  dynamic  loads  and  Ls  is  the  gear  mesh  life  based 
on  the  expanded  life  theory  and  static  loads  (as  illustrated  in 
fig.  4).  Comparing  figures  12  and  13  shows  that  the  gear  mesh 
life  decreased  when  the  maximum  dynamic  load  during  contact 
increased.  This  was  true  even  though  the  analysis  considered 


Figure  12.— Effect  of  speed  on  maximum  gear  tooth  dynamic  load  during 
contact  for  gear  data  from  table  I. 


the  load  along  the  complete  contact  length.  The  mesh  life  as 
a  function  of  speed  was  lowest  at  resonance. 

Effect  of  Mass,  Stiffness,  and  Damping  on  Gear  Life 

The  mass,  stiffness,  and  damping  of  a  gear  mesh  system 
significantly  affected  dynamic  load  and  life.  Modified 
TELSGE  was  run  using  the  mesh  data  in  table  I  while  varying 
the  equivalent  mass  A/^  and  keeping  all  other  parameters  the 
same.  The  life-speed  results  were  identical  when  plotted  on 
dimensionless  coordinates  (as  in  fig.  13).  Modified  TELSGE 
was  also  run  while  varying  the  equivalent  stiffness  and 
keeping  all  other  parameters  the  same.  Again,  the  life-speed 
results  were  identical  when  plotted  on  dimensionless 
coordinates  (as  in  fig.  13).  Thus  the  value  of  the  equivalent 
mass  or  the  equivalent  stiffness  had  no  effect  on  the  life-speed 
results  when  plotted  on  dimensionless  coordinates.  However, 
as  expected  from  equation  (31),  different  values  of  the 
equivalent  mass  or  the  equivalent  stiffness  produced  different 
values  for  the  resonant  speed.  The  equivalent  mass  and 
equivalent  stiffness  must  accurately  portray  the  gear  mesh 
being  modeled  for  the  calculated  resonant  speed  to  be  accurate. 

The  damping  force  in  the  dynamic  load  model  depends  on 
the  gear  system’s  viscous  friction  and  is  usually  an  unknown. 
Damping  ratios  f  between  0. 1  (in  eq.  (23))  and  0.2  were  used 
in  reference  11  to  correlate  analytical  and  experimental 
dynamic  load  gear  tests.  Here  damping  ratios  of  0. 10,  0.17, 
and  0.25  were  used  (fig.  14).  Decreasing  the  damping  ratio 
increased  the  dynamic  load  and  thus  shortened  the  mesh  life 
at  speeds  near  the  resonant  speed  and  one-half  the  resonant 
speed  (u/un  =  1.0  and  0.5,  respectively).  A  damping  ratio  of 
0.17  was  used  in  the  original  version  of  TELSGE  and  was 
used  in  this  study  for  all  other  figures. 

Effect  of  Speed  and  Contact  Ratio  on  Gear  Life 

Modified  TELSGE  was  used  to  predict  how  speed  and 
contact  ratio  affect  dynamic  load  and  gear  life.  Number  of 
teeth,  gear  size,  diametral  pitch,  pressure  angle,  and  gear  ratio 
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Dynamic  life  factor, 


TABLE  II,— DRIVER  GEAR  DATA 

[Set  I  used  for  baseline;  shaded  area  indicates  parameter  varied  from 
baseline.] 


•17V\ 

^°\V\ 


Set  Number  Pitch  Diametral  Pressure  Gear  Contact 

of  radius,  pitch  angle,  ratio  ratio 

teeth  cm  deg 


Figure  16.— Effect  of  contact  ratio  on  gear  mesh  life  for  parametric  study 
data  from  table  fl. 


For  a  mesh  with  a  contact  ratio  of  1 .6,  two  pairs  of  teeth  are 
in  contact  60  percent  of  the  time  and  one  pair  is  in  contact 
40  percent  of  the  time.  Low-contact-ratio  gears  have  contact 
ratios  between  1  and  2.  In  the  parametric  studies  the  contact 
ratio  ranged  from  1.41  to  1.92. 

For  the  data  in  figure  15  the  dynamic  life  factor  was  plotted 
as  a  function  of  contact  ratio  in  figure  16  for  speeds  w/wn  of 
0.7  and  1.0.  A  sixth-order  polynomial  curve-fit  was  used  to 


generate  the  curves.  At  a  constant  normalized  speed  the 
dynamic  life  factors  were  about  the  same  for  meshes  with 
contact  ratios  between  1 .4  and  1 .6  but  were  significantly  higher 
for  meshes  with  higher  contact  ratios. 

With  higher  contact  ratios  the  equivalent  stiffness  (fig.  9) 
had  a  smaller  duration  of  single-tooth-pair  contact  and  thus 
a  smoother  transition  of  double-  to  single-  to  double-tooth- 
pair  contact.  This  resulted  in  lower  dynamic  load  factors  and 
higher  dynamic  life  factors.  For  the  sets  studied,  the  resonant 
speed  varied  with  equivalent  mass,  mean  equivalent  stiffness, 
pressure  angle,  and  number  of  teeth. 

A  general  design  chart  for  the  dynamic  life  factor  of  a  gear 
mesh  was  developed  from  the  parametric  studies  (fig.  17).  The 
objective  was  to  determine  the  dynamic  life  factor  as  a  single 
simple  function  of  speed  and  contact  ratio  to  be  used  when 
a  computer  and  program  TELSGE  are  not  available.  The  heavy 
solid  line  represents  the  best  fit  of  the  results  of  the  parametric 
studies.  For  co/cun  <  0.5  the  dynamic  life  factor  can  be  read 
directly  from  the  plot  by  using  the  scale  on  the  left.  For 
<o/a>„  >  0.5  the  dynamic  life  factor  is  the  product  of  the  value 
of  the  curve  (using  the  scale  on  the  right)  and  the  contact  ratio 
to  the  sixth  power.  The  light  dotted  lines  represent  the  actual 
results  of  the  parametric  studies  and  indicate  the  possible  error 
when  using  the  chart.  An  example  problem  given  in  appendix 
B  demonstrates  the  use  of  the  design  chart.  A  simplified  hand 
calculation  of  gear  tooth  stiffness  is  also  given  in  appendix  B. 


Figure  17.— Dynamic  life  factor. 
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Summary  of  Results 

How  dynamic  load  affects  the  pitting  fatigue  life  of  external 
spur  gears  was  predicted  by  using  a  modified  version  of  the 
NASA  computer  program  TELSGE  to  perform  parametric 
studies.  TELSGE  was  modified  to  include  a  surface  pitting 
fatigue  life  analysis.  The  parametric  studies  modeled  low* 
contact-ratio  involute  gears  with  no  tooth  spacing  or  profile 
errors.  The  following  results  were  obtained: 

1 .  Gear  life  predictions  based  on  dynamic  loads  generally 
differed  significantly  from  those  based  on  static  loads  and  were 
strongly  influenced  by  the  maximum  dynamic  load  during 
contact. 

2.  Gear  mesh  operating  speeds  strongly  affected  predicted 
dynamic  loads  and  thus  gear  life.  In  most  cases  studied,  meshes 
operating  at  a  resonant  speed  or  one-half  the  resonant  speed 
had  significantly  shorter  lives  than  the  life  based  on  static  loads. 
Meshes  operating  above  resonance  had  significant  longer  lives. 


3.  In  general,  meshes  with  higher  contact  ratios  had  higher 
predicted  dynamic  life  factors  than  meshes  with  lower  contact 
ratios. 

4.  Damping  significantly  affected  predicted  gear  mesh  life 
for  meshes  operating  at  or  near  a  resonant  speed  or  one-half 
the  resonant  speed. 

5.  A  solution  for  dynamic  load  converged  with  only  a  few 
iterations  of  gear  tooth  stiffness  and  dynamic  load  because  the 
Hertzian  deflection  was  relatively  small  in  comparison  with 
the  total  gear  tooth  deflection. 

Lewis  Research  Center 

National  Aeronautics  and  Space  Administration 
Cleveland,  Ohio,  April  3,  1986 
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Appendix  A 
Symbols 


B  material  constant  (2.23  x  108  N/m 1 979,  35  000  lb/in. >m, 
ref.  21) 

Qq  equivalent  damping  per  unit  face  width,  N  sec/m2 
(lb  sec/in.2) 

Cv  dynamic  life  factor 

c  contact  ratio 

d  distance  of  inscribed  parabola,  m  (in.) 

E  modulus  of  elasticity,  Pa  (psi) 

e  Weibull  exponent 

/  gear  tooth  face  width,  m  (in.) 

ho  tooth  thickness  at  root  radius,  m  (in.) 

/o  beam  cross-sectional  moment  of  inertia,  m4  (in.4) 

/  contact  position  index 

J  number  of  intervals;  or  polar  mass  moment  of  inertia 
per  unit  face  width  (1/2  mrb  for  disk),  kg  m  (lb  sec2) 

j  interval  index 

K  combined  stiffness  per  unit  face  width.  Pa  (psi) 

K eq  equivalent  stiffness  per  unit  face  width.  Pa  (psi) 

K eq  mean  equivalent  stiffness  per  unit  face  width,  Pa  (psi) 
k  gear  tooth  stiffness  per  unit  face  width.  Pa  (psi) 

L  life  for  90-percent  probability  of  survival,  Mrev 

f  involute  length,  m  (in.) 

M  effective  mass  per  unit  face  width,  7/r^,  kg/m 
(lb  sec2/in.2) 

equivalent  mass  per  unit  face  width,  kg/m  (lb  sec2/in.2) 
m  mass  per  unit  face  width,  kg/m  (lb  sec2/in.2) 
ma  module,  mm/tooth 

N  number  of  teeth 

P  diametral  pitch,  teeth/in. 

Pd  dynamic  load  per  unit  face  width,  N/m  (lb/in.) 

Ps  static  load  per  unit  face  width,  N/m  (lb/in.) 
pb  base  pitch,  2 *rb/N,  m  (in.) 

Q  normal  load,  N  (lb) 

Q  average  interval  load,  N  (lb) 

Q,  tangential  load,  Q  cos  <0,  N  (lb) 

R  radius  of  curvature,  m  (in.) 
rb  base  radius,  rp  cos  <p,  m  (in.) 

r0  outside  radius,  m  (in.) 

rp  pitch  radius,  Nm0/2  =  N/2P,  m  (in.) 


rr  root  radius,  m  (in.) 

s  displacement,  m  (in.) 

T  torque,  N  m  (lb  in.) 

t  tooth  thickness  at  pitch  radius,  m  (in.) 

X  relative  displacement,  m  (in.) 

jc  contact  position,  m  (in.) 

Ax  interval  length,  m  (in.) 

Y  Lewis  form  factor 

Z  contact  length,  m  (in.) 

Z\  contact  length  from  pitch  point  to  start  of  contact,  m  (in.) 

Z2  contact  length  from  pitch  point  to  end  of  contact,  m  (in.) 

a  pressure  angle  at  root  radius,  deg 

y  density,  kg/m3  (lb/in.3) 

d  beam  deflection,  m  (in.) 

f  clamping  ratio 

7}  life  for  90-percent  probability  of  survival,  millions  of 

stress  cycles 

6  angular  displacement,  rad 

Ep  curvature  sum,  m_I  (in.^1) 

£p  average  curvature  sum,  m~l  (in.-1) 

<p  pressure  angle,  deg 

o)  speed,  rpm 

resonant  speed,  rpm 

Subscripts: 
d  dynamic  life 

/  gear  tooth  stiffness-dynamic  load  iteration  index 
/  contact  position  index 

j  interval  index 

m  mesh 

max  maximum  during  contact  position 
s  static  life 

t  tooth 

1  driver  gear 

2  driven  gear 

Superscripts: 

(I)  first  pair  of  teeth  in  contact 

(II)  second  pair  of  teeth  in  contact 
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Appendix  B 
Example  Problem 


Example  problem.— Dt termine  the  dynamic  life  factor  of 
the  mesh  from  the  data  given  in  table  III  when  the  driver  gear 
is  rotating  at  5000  rpm. 

Solution.— The  pitch  radii  are 


Pt  = 


2xrb  i  _  2tt(6.138  cm) 


Nx  32  teeth 

From  equation  (33)  the  contact  ratio  is 


=  1.205  cm  (0.475  in.) 


1  1  /  ,  1  cm  \ 

r_  i  —  -  =  -  (32  teeth)(  4.233  mm/tooth  x - ] 

pl  2  2  \  10  mm/ 


=  6.773  cm  (2.667  in.) 


(4! 


r n  y  =  -  (100  teeth)!  4.233  mm/tooth  x 


1  cm 


10  mm  j 


=  21  165  cm  (8  333  in.) 


The  base  radii  are 

rh  i  =rp  ,  cos  <  =  (6.773  cm)  cos  25°  =6. 138  cm  (2.417  in.) 

rh  :  =  (21.165  cm)  cos  25*  =  19.182  cm  (7.552  in.) 

From  equations  (5)  and  (6)  the  contact  lengths  from  the  pitch 
point  to  the  start  and  end  of  contact  are 


1.858  cm  _ 

/*  — _ _ — i  —  1  vl 

1.205  cm 


The  masses  per  unit  face  width  of  the  driver  and  driven  gears 
can  be  approximated  by 


m,  =  yxTpj  =  (7833  kg/m3)* 


=  112.886  kg/m  (1.637  x  10“2  lb  sec2/in.2) 


6.773  cm  x 


1  m 
100  cm 


m2  =  (7833  kg/m 


21.165  cm  x 


1  m  V 
100  cm/ 


=  1102.337  kg/m  (1.598  X  HP1  lb  sec2/in.2) 


C|  =  \  (2 1 .590  cm)2  -  ( 19. 182  cm)2 

-  (21.165  cm)  sin  25°=  0.964  cm  (0.379  in.) 
r:  =  V(7. 196  cm)2  -  (6.138  cm)2 

-  (6.773  cm)  sin  25°=  0.894  cm  (0.351  in.) 

From  equation  (7)  the  contact  length  is 
Z  =  (0.964  cm)  +  (0.894  cm)  =  1.858  cm  (0.730  in.) 
The  base  pitch  is 


The  effective  masses  per  unit  face  width  are 

l  2 

-mrti 

M,  =  -r-  = - =  -m,  =  -(112.886  kg/m) 

4.1  dA  2  2 

=  56.443  kg/m  (8.185  x  10~3  lb  sec2/in.2) 

M2  =  ^(1 102.337  kg/m)  =  551.169  kg/m 
(7.990  x  10  ”2  lb  sec2/in.2) 


TABLE  III. -GEAR  MESH  DATA  USED  IN  DYNAMIC  LIFE 
FACTOR  EXAMPLE  PROBLEM 


Parameter 

Driver  gear 

Driven  gear 

Number  of  teeth 

32 

100 

Outside  radius,  cm  (in.) 

7.196  (2.833) 

21.590  (8.500) 

Root  radius,  cm  (in.) 

6.246  (2.459) 

20.638  (8.125) 

Lewis  form  factor 

0.433 

0.521 

Module,  mm/tooth  (Pitch,  teeth/in.) 
Face  width,  cm  (in.) 

Pressure  angle,  deg  (rad) 

Tooth  thickness  at  pitch  radius, 
cm  (in.) 

Modulus  of  elasticity.  Pa  (psi) 

4.233  (6) 
6.350  (2.500) 
25  (0.436) 
0.665  (0.262) 

2.068  x  10n  (30  x  I06) 

Density,  kg/m3  (lb/in.3) 

7833  (0.283) 

From  equation  (22)  the  equivalent  mass  per  unit  face  width  is 

(56.443  kg/m)(55 1.169  kg/m) 

”  (56.443  kg/m)  +  (551.169  kg/m) 

=  51.200  kg/m  (7.424  x  10  ~3  lb  sec2/in.2) 

Determining  teeth  stiffnesses  by  the  methods  of  Cornell 
(ref.  17)  requires  the  use  of  a  computer.  For  this  example  the 
stiffness  calculations  will  be  simplified  by  modeling  the  gear 
teeth  as  cantilever  beams  of  uniform  strength  (beams  in  which 
the  section  modulus  varies  along  the  beam  in  the  same 
proportion  as  the  bending  moment).  The  pressure  angles  at 
the  root  radii  of  the  driver  and  driven  gears  (fig.  IE)  are 


Figure  18  — Gear  tooth  model  for  simplified  stiffness  calculations. 


a,  =  cos 


-1  frb.\\  -I  ^6  138  cm\ 

1  (  —  I  =  cos  I - ) 

\rrl/  \6.246  cm/ 


=  10.670°  =0.186  rad 


<*2  =  COS 


/l9.]82cin\  =  21.651°  =  0.378  rad 
\20.638  cm/ 


From  reference  22  the  tooth  thickness  at  the  pitch  radius  is 
related  to  the  tooth  thickness  at  the  root  radius  by 


=  2r  r** 

V r. 


+  (tan  a  -  a)  -  (tan 


<S>  ~  *>)j 


Therefore  the  teeth  thicknesses  at  the  root  radii  are 

V i  -  2 rr  l  — ^ - (tan  a,  -  ot,)  +  (tan  <p  -  <p) 

L2rp.i 

[  0.665  cm 

- (tan  0.186  -  0.186) 

2(6.773  cm) 

+  (tan  0.436  -  0.436)J 


[  0.665  cm 

- (tan  0.378-0.378) 

2(21.165  cm) 


+  (tan  0 


.436  -  0.436) 


=  1.095  cm  (0.431  in.) 


From  reference  23  and  figure  18  the  distance  of  the  inscribed 
parabola  is 


where  Y  is  the  Lewis  form  factor.  Thus  the  distances  of  the 
inscribed  parabolas  are 


,2  (0.960  cm)2 

ftp.i  _ 

6 m0Y\  /  1  cm  \ 

6(  4.233  mm/tooth  x - 1(0.433) 

V  10  mm/ 


=  0.838  cm  (0.330  in.) 


(1.095  cm)2 


6(  4.233  mm/tooth  x-?-^-)(0.521) 
V  10  nun/ 


=  0.906  cm  (0.357  in.) 


The  inscribed  parabola  in  figure  18  is  a  cantilever  beam  of 
uniform  strength.  From  reference  24  the  deflection  for  the 
beam  is 


6  =  where  l°  = 

3£7q  12 

The  gear  tooth  stiffness  per  unit  face  width  is 

Q,  _  3 EIp  =  Ehp 
6f  2 dlf 

For  the  driver  and  driven  gears »  respectively. 


(2.068  x  10 


ii  /  1  m  V 

11  Pa)(  0.960  cm  x - 

\  100  cm/ 

/  1  m  V 

(  0.838  cm  x - 

\  100  cm/ 


=  0.960  cm  (0.378  in.) 


=  3.886  x  10 10  Pa  (5.636  x  106  psi) 
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/  lm  V 

8(  0.916  cm  x  — - ) 

\  100  cm/ 


=  4.564  x  10 10  Pa  (6.599  x  106  psi) 


0  Z  -  Pu  Pb  z 

Contact  position,  x 

Figure  19.— Equivalent  gear  tooth  stiffness  as  a  function  of  contact  position. 


From  equation  (25)  the  combined  stiffness  per  unit  face  width  From  equation  (31)  the  resonant  speed  in  terms  of  driver  gear 


rotations  is 


(3.886  x  IQ10  Pa)(4.564  x  IQ10  Pa) 

~~  (3.886  x  10l°  Pa)  +  (4.564  x  1010  Pa) 

=  2.099  x  1010  Pa  (3.040  x  106  psi) 


For  this  example  it  is  assumed  that  the  combined  stiffness 
is  constant  with  respect  to  contact  position.  During  single¬ 
tooth-pair  contact  the  equivalent  stiffness  is  K .  During  double¬ 
tooth-pair  contact  the  equivalent  stiffness  is  2 K.  The  mean 
equivalent  stiffness  per  unit  face  width  (fig.  19)  is  given  as 


*eq  =  -|  K^dX 
0 


i/n 


K^dx+  K^dx+  K^dx 


=  i(2#r[(Z  “  p»)  ~  °]  +  K\pb  -  (Z  -  pb)\ 

+  2K[Z-Pb]} 

=43-2(z)]‘K5^) 

=  (2.099  x  10'°  Pa)  ^3  - 
=  3.571  x  I0'°  Pa  (5.172  x  106  psi) 


V— COS  <p 

__Meq__/60\ 
(V,  W 


'3.571  x  10'°  Pa 
51.200  kg/m 

32  teeth 


cos  25“ 


©■ 


7143  rpm 


At  a  driver  operating  speed  of  5000  rpm 


a)  5000  rpm 
o)n  7143  rpm 

From  figure  17  for  =  0.70 


=  0.70 


£ = 004 

c6 


and  the  dynamic  life  factor  is 

Cv  =  0.04  c6  =  0.04  (1.54)6  =  0.53 

Thus  about  a  50-percent  decrease  in  life  compared  with  that 
using  static  loads  is  predicted  for  this  example.  Note  that  the 
simplified  stiffness  model  used  in  the  example  may  produce 
erroneous  values  for  the  resonant  speed.  The  mean  value  of 
the  equivalent  stiffness  per  unit  face  width  for  this  example  was 
computed,  by  using  Cornell’s  method  (ref.  17)  and  TELSGE, 
as  2.741  x  10 10  Pa  (3.975  x  106  psi).  This  produced  a 
resonant  speed  of  6260  rpm  and  a  dynamic  life  factor  of  0.56. 
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